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Abstract. The Coma cluster is one of the richest known 
cluster of galaxies, spanning about 4 dex in density. Hence 
it is the ideal place to study the structure of galaxies as 
a function of environmental density in order to constrain 
the theories of galaxy formation and evolution. For a mag- 
nitude limited sample of 35 E and SO galaxies we obtained 
long slit spectra to derive the rotation curves, the velocity 
dispersion profiles and the radial gradients of the Mg, Fe 
and H/3 line indices. Here we report on two early - type 
galaxies which turned out to host the largest kinematically 
peculiar cores yet found in "normal" early - type galax- 
ies: NGC 4816 hosts a decoupled counter rotating core 
with a radial extension along the major axis of 2.7 kpc, 
while IC 4051 has a co- rotating peculiar core with a sizes 
of 3.4 kpc. We combine our data with HST photometry 
and show that both cores are flattened central stellar disks 
which contribute less than 1 % to the total V band light 
of the galaxies, but are nevertheless conspicuous (1-2 
x 10 9 L ). The metallicity of the cores is 0.25 dex super 
solar and drops to solar and sub solar in the outer part 
of NGC 4816 and IC 4051, respectively. The mean stellar 
population in both central disks is old (8 - 12 Gyr) and 
highly overabundant in Mgb relative to <Fe> (« 0.5 dex). 
We discuss the evidence that these central disks formed 
via dissipational major merger events. 



Key words: galaxies: elliptical and lenticular, cD - kine- 
matics and dynamics - abundances - formations - indi- 
vidual: NGC 4816 - IC4051 



Susm.um- 



Send offprint requests to: D. Mehlert (email: mehlert<l 
muenchen.de) 

* Visiting astronomer at the German- Spanish Astronomi- 
cal Center, Calar Alto, operated by the Max-Planck-Institut 
fur Astronomie, Heidelberg jointly with the Spanish National 
Commission for Astronomy. 

** Visiting astronomer at the Michigan - Dartmouth - M.I.T. 
Observatory, Kitt Peak, Arizona, operated by a consortium of 
the University of Michigan, Dartmouth College and the Mas- 
sachusetts Institute of Technology. 



1. Introduction 

Kinematically peculiar cores are generally understood as 
being fossil fingerprints of the merging - and thus for- 
mation - history of the host galaxies (Kormendy 1984, 
Franx & Illingworth 1988, Bender 1988, Jedrzejewski & 
Schechter 1988). They have mostly been explained by ei- 
ther accretion of small spheroidal satellites by larger el- 
lipticals (Kormendy 1984, Balcells & Quinn 1990, Balcells 
1991) or as the result of major merging events by star 
dominated systems (Schweizer 1990, Barnes 1992a, Ben- 
der & Surma 1992, Hernquist 1993, Heyl et al. 1994). Up 
to now kinematically decoupled cores with radial exten- 
sions between 0.2 - 2.7 kpc have been found in a number 
of elliptical galaxies - lying in the field and in groups as 
well as in clusters (Bender 1990; see Sect. 5 for a detailed 
list). 

Most peculiar cores have been shown to be due to a 
rapidly spinning disk - like component which presumably 
formed dissipatively (Bender 1988, 1990, Franx & Illing- 
worth 1988, Rix & White 1992, Surma & Bender 1995). 
As shown by Bender & Surma (1992) the high absorp- 
tion line strengths of the decoupled regions further argue 
against accretion of compact small galaxies and suggest a 
scenario analogous to massive spiral - spiral mergers (see 
also Davics ct al. 1993). In this model the prc-enriched gas 
in the two merging components lost its angular momen- 
tum, settled in the central region and underwent a strong 
star formation phase. Though spiral - spiral mergers may 
still form ellipticals today (Kormendy & Sanders 1992), 
it is unlikely that most decoupled cores formed recently. 
Their location in the Mg - a diagram indicates that the 
corresponding ellipticals have ages similar to normal el- 
lipticals. Recent HST images of kinematically decoupled 
cores (Carollo et al. 1997) show that 45 - 60 % of them 
are in fact stellar central disks and confirm the idea, that 
most of them indeed have formed by dissipational merging 
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processes. 

Finally HST imaging by Jaffe et al. (1994), Lauer et al. 
(1995) & Faber et al. (1997) has shown that the most lu- 
minous, anisotropic, boxy galaxies possess surface profiles 
with shallow cores while low luminous, rotationally flat- 
tened, disky ellipticals have power-law profiles (confirming 
earlier claims by Nieto et al. 1991). 

Here we present kinematical and line index profiles as 
well as HST surface brightness profiles for the two Coma 
galaxies IC 4051 and NGC 4816, which host the largest 
cores yet found in early - type galaxies. 
In Sect. 2 we describe the spectroscopic observations, data 
reductions and the results of the kinematical analysis and 
line strength measurements. The photometric data ob- 
tained with the Hubble Space Telescope (HST) and the 
MDM 1.3 m telescope are described in Sect. 3, where we 
derive the surface brightness profiles and other photomet- 
ric parameters. We summarize the combined spectroscopic 
and photometric results in Sect. 4. Finally we discuss for- 
mation scenarios for NGC 4816 and IC 4051 in Sect. 5. 

2. Spectroscopy 

2.1. Observations & data analysis 

Long-slit spectra centered on the 5170 A Mg triplet 
of NGC 4816 and IC 4051 were taken along their major 
axes. NGC 4816 was observed at the 3.5 m telescope of the 
German-Spanish Astronomical Center on Calar Alto (CA) 
in 1996 (see Table @). IC 4051 was observed with the 2.4 m 
telescope of the Michigan-Dartmouth-M.I.T. (MDM) ob- 
servatory on Kitt Peak in 1995. Additionally, during both 
runs we obtained spectra for ten K - Giant template stars, 
trailed and wiggled across the slit. The standard reduction 
(bias & dark subtraction, flat fielding, wavelength calibra- 
tion, sky subtraction, correction for CCD misalignment) 
was carried out with the image processing package MI- 
DAS provided by ESO. 

After continuum subtraction and radial binning along 
the slit we achieved a S/N > 40/ A for all of our data. 
We determined the line-of-sight-velocity-distributions 
(LOSVDs) by using the Fourier Correlation Quotient 
(FCQ) method (Bender 1990), which provides us with the 
stellar rotational velocities v rot , the velocity dispersions 
cr and the first orders of asymmetric (H 3 ) and symmetric 
(H4) deviations of the LOSVDs from real Gaussian pro- 
files (van der Marel & Franx 1993 ; Bender et al. 1994). As 
expected, we find that the FCQ method is little influenced 
by template mismatching (Bender 1990). The "best" tem- 
plate star minimizing the asymmetric deviations (H3) of 
the LOSDV from a real Gaussian profiles turned out to 
be HR 6817. Following Bender et al. (1994), Monte Carlo 
simulations were made to find the best fit-order for the 
continuum, to check for systematic effects and to estimate 
errors. For NGC 4816 two exposures with an atmospheri- 



cal seeing of FWHM = 1.2" and one with FWHM = 3.7" 
were available. The latter one was excluded for the final 
analysis, but was used to check, whether the large slit 
width (3.6") affected the spectral resolution significantly 
for good seeing conditions. The three individual exposures 
did not produce significant differences within the errors. 
Line strength indices were derived following Faber et al. 
(1985) and Worthey (1992) from flux calibrated spectra, 
rebinned radially as before. Furthermore, we corrected for 
velocity dispersion broadening and calibrated our mea- 
surements to the Lick system using stars from Faber et al. 
(1985). The errors are derived from photon statistics and 
CCD read-out noise. 
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Fig. 1. (a) The measured Gauss width a at each wavelength 
for one single frame of NGC 4816. We measured a inside a 
radius of 7", which turned out to be the radius at which the 
focus variation starts to effect the continuum of the measured 
spectra. For H^, Mg;,, Fes and Fe6 the line windows defined 
by Faber et al. (1985) are indicated (solid lines). For the Mgi 
and Mg2 indices the line windows (long dashed lines) and the 
windows for their pseudo continuum (short dashed lines) are 
shown, (b) The spectra of NGC 4816 at different radii before 
(thin solid line) and after (fat dashed dotted line) the focus 
correction. The different continua shapes in the uncorrected 
spectra are evident. Like in Fig. (a) the windows for the mea- 
sured indices are indicated. 

The line index profiles for NGC 4816 needed further treat- 
ment. At small radii some of them show unreal features 
caused by the varying focus of the spectrograph in the 
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Table 1. Setup of the observing runs 



dispersion direction, which depends slightly on the wave- 
length. Following Gonzales (1993), this effect is only de- 
tectable if the variation of the focus' point spread function 
(PSF) is dominant compared to the atmospherical seeing 
(FWHM). For NGC 4816 the focus variation was the dom- 
inant effect and we developed a new method to correct for 
this "focus effect": For each single exposure of NGC 4816 
we measured the central Gaussian width cr(A) of the ra- 
dial profile at each wavelength (Fig. 0a) and noted that 
er is larger at the blue and red end of the spectra com- 
pared to the central wavelength. Therefore blue and red 
light is spread more radially than light in the center of the 
spectrum. Consequently, the spectrum's continuum shape 
varies with radius (see Fig. thin solid line) and leads to 
unreal central features in the measured line index profiles. 
The Mgi and Mg 2 indices defined by Faber et al. (1985) 
are affected most, because their pseudo continua are rj 
200 A apart from the line windows. For all the other line 
indices the continuum and line windows are close to each 
other. and the Fe indices are in the red and blue wing of 
the spectra and thus slightly affected, while the Mgf, index 
at the central wavelength is almost unaffected. To correct 
for this focus variation in the exposures of NGC 4816 we 
convolved all radial profiles with the wavelength depend- 
ing width a broad = (<7? nax - a 2 {\)) 1/2 (see Fig. |l|a). As a 
result the radial profiles have the width a ma x = 2.47" at 
all wavelength. Furthermore, the continuum shapes of the 
spectra are the same for all radii (see Fig. [|b; fat dashed 
dotted line) and the artificial features in the line index 
gradients disappear (see Mehlert 1998 for a more detailed 
description). For the spectra of NGC 4816 the maximum 
broadening width is o%™ d = (a 2 max - <J 2 mm ) 1/2 =1-96" 
{o~min =1.5" ; see Fig. [lja). The resulting PSF for the mea- 
sured line index profiles (shown in Fig. ||b) has <r res = 

((^f) 2 + Seeing) 112 = 2 -°" ^ FWHM res = 4.7" . 

For IC 4051 no "focus correction" was necessary, because 
the atmospherical seeing was dominant during the obser- 
vations. 



2.2. Kinematical profiles and line index gradients 

Figs. ^ and || show the kinematical profiles (a) together 
with the line indices (b) of NGC 4816 and IC 4051, re- 



spectively. 
NGC 4816: 

This galaxy contains a kinematically decoupled core, 
which is counter-rotating with respect to the main 
spheroid at a maximum velocity v ma x = 30 km/s. The 
kinematical break radius is rj, = 4", which is 2.7 kpc in 
the Coma cluster (1' = 40 kpc; H Q = 50 km/s/Mpc). The 
velocity dispersion is o\, = 220 km/s at r;, and its gradi- 
ent gets much steeper inside the core, rising to a central 
value of a a « 260 km/s. Finally, H3 changes sign at the 
core radius. The central velocity dispersion measured by 
Davies et al. (1987) lies a bit above our profile, but still 
agrees with it within the errors. 

For NGC 4816 the line indices measured with (large 
symbols) and without (small dots) the focus correction 
described in Sect. 2 are shown and the need and improve- 
ment for the Mg2 index is obvious. Without focus correc- 
tion the profile of the Mg2 index is disturbed out to 7" . 
For the focus corrected profiles we find an increase of slope 
in the Mg'2 as well as in the Mg& line profile inside 8", 
which is larger than if, (see section |]for explanation). Be- 
yond 8" the Mgh index is constant. The mean iron in- 
dex <Fe> = (Fe5 + Fee) / 2 and the combined index 
[MgFe] = (Mg h <Fe>) 1 / 2 do also increase inside 8" and 
stay constant beyond that point. The index shows a 
slight increase towards the outer parts of the galaxy, but is 
dropping at the very outer (less reliable) data point. The 
Mg2 value derived by Davies et al. (1987) is a little larger 
than our data, but still agrees within the errors. 
IC 4051: 

This galaxy contains a peculiar core which is co-rotating 
with the main spheroid at a maximum projected velocity 
Vmax — 60 km/s. The core's kinematics are detectable to 
a radius of 5 ", which is 3.4 kpc in the Coma cluster. It 
is not clear whether the core of IC 4051 is really decou- 
pled from the rotation of the main body, since the velocity 
curve does not drop to exactly at the break radius r^. 
The velocity dispersion is en, = 220 km/s at r?, and its gra- 
dient gets much steeper inside the core, rising to a central 
value of cr « 280 km/s. Note that the H3 parameter is 
negative for the co-rotating core. J0rgensen et al. (1995) 
as well as Davies et al. (1987) measured the central ve- 
locity dispersion for this galaxy, deriving values in good 
agreement with ours. 



4 



Mehlert et al.: The kinematically peculiar cores of NGC 4816 and IC 4051 



NGC 4816; r e = 55.75" ; PA = 78° (MA) 
a) b) 




10 20 30 10 20 
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Fig. 2. (a) Kinematical parameters and measured line indices (b) for NGC 4816 along the major 
axis. The curves are folded around the nucleus of the galaxy and different symbols refer to different 
sides. For (b) the profiles before (small dots) and after focus correction (large symbols) are shown. 
Open squares are results from Davies et al. (1987) who used a 4" x4" aperture. We determined the 
equivalent radius along our used slit via r d87 =0.5 x (4"x4"/ slitwidth). Panels from bottom to 
top: (a) (1) Rotation velocity v(r) relative to the system velocity, (2) velocity dispersion er(r), (3-4) 
Gauss- Hermite moments H3 and H4. (b) (1-2) Mg2 and Mgt, (3) mean Fe index <Fe> = (Fes + 
Fe 6 )/2, (4) the combined index [MgFej = (Mg 6 <Fe>) 1/2 , (5) R fJ index. 



The Mg;, index is unusually high, showing values up to 
5.5 A in the center of the galaxy. It is constant in the core 
region, but decreases strongly beyond r^. The asymmetry 
of the Mgh profile in the center of the galaxy as well as 
in the outer parts is probably due to contamination from 
[NI] emission lines in the interstellar medium (Goudfrooij 
& Emsellem 1996). [NI] emission can artificially increase 
the measured Mgb index. The asymmetry itself arises from 
the fact that the stars and the gas may have different ve- 
locities up to 100 km/s (Bcrtola et al. 1995). Additionally, 
as shown by Goudfrooij & Emsellem (1996), the Mg 2 in- 



dex is much less disturbed by the [NI] emission lines, be- 
cause of its different continuum definition. However, the 
S/N and spectral resolution of our data are insufficient to 
quantify this effect. In this galaxy the mean Iron index 
<Fe> shows a steep increase towards the center starting 
at the break radius, while the combined index [MgFe] is 
constantly increasing towards the center of the galaxy. A 
steep decrease of the absorption index is present out- 
side 7" . The inspection of the absorption lines shows 
that they might be partially filled by emission like features 
outside this radius. As above, the S/N and resolution of 
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IC 4051: r =15.56' 
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Fig. 3. Kinematical parameters (a) and measured line indices (b) for IC 4051 along the major axis. 
The curves are folded around the nucleus of the galaxy with different symbols referring to different 
sides. Open triangles represent data from j0rgensen et al. (1995) and open squares those from Davies 
et al. (1987). The latter used a 4" x4" aperture. We determined the equivalent radius along our used 
slit via r<287 = 0.5 x (4" x4"/ slitwidth). Panels from bottom to top: (a) (1) Rotation velocity v(r) 
relative to the system velocity, (2) velocity dispersion <r(r), (3-4) Gauss-Hermite moments H3 and 
H 4 . (b) (1-2) Mg 2 and Mg b , (3) mean Fe index <Fe> = (Fe 5 + Fe 6 )/2, (4) the combined index 
[MgFe] = (Mg b <Fe>) 1/2 , Hp index. 



the data do not allow to quantify the effect. Central values 
of Mg2 are provided again by J0rgensen et al. (1995) and 
Davies et al. (1987) and agree with our data. 



3. Photometry 

3.1. HST observations & data analysis 

HST images were available for both galaxies, NGC 4816 
(Principal Investigator: John Lucey; Proposal ID: 5997) 
and IC 4051 (Principal Investigator: James Westphal; 



Proposal ID: 6283). For NGC 4816 exposures with 800 s 
in total are available with the F606W Filter, while for 
IC 4051 eight full-orbit WFPC2 exposures (20500 s in to- 
tal) with the F555W Filter (A c = 5957 A) and two 
(5200 s in total) with F814W (A c = 7940 A) were ob- 
tained. Cosmic ray elimination and the production of a 
median frame from multiple observations were carried out 
within " STSDAS" packages of the image processing sys- 
tem IRAF § (Tody 1993). 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 



(i 



Mehlert et al.: The kinematically peculiar cores of NGC 4816 and IC 4051 



The subsequent analysis was done within MIDAS using 
the data of the Planetary Camera (PC), which has a scal- 
ing of 0.04555" /pixel. The sky values were determined 
from the Wide Field Camera (WFC) images in areas which 
were at least 55" away from the galaxy's center. After sky 
subtraction the foreground stars and galaxies were either 
replaced with values from the opposite side of the galaxy 
or totally removed and excluded from the analysis. The 
photometric zero points were derived and converted to 
the Johnson-Cousins system following the procedure de- 
scribed in Holtzman ct al. (1995). We added a constant 
shift of 0.115 mag to correct for infinite aperture and de- 
termined the I(F814W), V(F606W)and V(F555W) band 
zero points to be 21.67 mag, 23.33 mag and 22.57, respec- 
tively. 

The deviations of the isophotes from ellipses were 
studied by means of Fourier series expansions (Bender & 
Mollenhoff, 1987). The resulting fourth cosine coefficients 
<i4 indicates boxy (negative values) or disky isophotes 
(positive values). Additionally, we estimated the bright- 
ness profile of the existing central disks with the decom- 
position method developed by Scorza & Bender (1995). 
The principal idea is to consider galaxies with a positive 
d4 parameter to be a composition of an ideal elliptical 
bulge with 04 = and a thin disk. The model disk is de- 
scribed by an exponential intensity profile and parameter- 
ized by the scale length xu, the central surface brightness 
Jo, the position angle P.A. and the inclination anglecl i. 
These free parameters are varied in an iterative process 
until a bulge having perfectly elliptical isophotes (04 « 0) 
remains. We also determined the ratio of the disk's and 
bulge's luminosity (D/B = ~Ld/(Lg - L.d)) for the whole 
galaxy and inside the kinematically detected break radius 
ib- For the distance module of the Coma cluster we used 
(m — M) = 35.7 mag and the V and I absolute magnitude 
of the Sun are M Q y = 4.84 mag and Mqj = 4.03 mag. 

For IC 4051 we also derived the V-I colour and its 
gradient. Since the PSF of the HST is different for the 
two filters and varies with the position on the CCD, we 
convolved each frame with the PSFEl of the other filter. 
Finally we get the same resulting PSF for both frames. 

The values of the total magnitude mj-, the surface 
brightness SB and the colour have been corrected for 
galactic extinction an if - corrected. We have used the 
B band galactic extinction determined by Burstcin et al. 
(1987) for both NGC 4816 (A B = 0.03 mag) and IC 4051 
(Ab — 0.05 mag). For the filters we used, we assumed A\ 
= k\As and following Seaton (1979) we determined the 
values k y = 0.79, k R = 0.58 and kj = 0.48. For the K- 
correction we used values provided by Rocca-Volmerange 
& Guiderdoni (1988) {K v = 0.041 mag, K R = 0.028 mag, 



versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation of the USA. 

2 where i — 0° means face on, i = 90° edge on. 

3 The PSFs where obtained with the TINYTIM routine. 



Ki = 0.0 mag). Finally, we corrected the surface bright- 
ness for cosmological dimming. 



3.2. MDM observations & data reduction 

We also obtained Kron-Cousins R band CCD surface pho- 
tometry for NGC 4816 at the 1.3 m telescope of the MDM 
observatory in March 1995. We used a Loral 2048 x 2048 
chips with a scaling of 0.508" per pixel and observed NGC 
4816 for 500s with a PSF of 2". Sky subtraction, cos- 
mic ray elimination and the removal of foreground stars 
were performed with MIDAS packages. We carried out the 
same isophote analysis and surface brightness fit as for 
the HST data. The zero point was estimated by using V 
band aperture photometry from Burstein et al.(1987) and 
using a typical colour for elliptical galaxies V-R = 0.65 
mag (Poulain & Nieto 1994). Total magnitudes hit and 
half luminosity radii r e were derived using the algorithm 
of Saglia et al. (1997a), which fits the surface brightness 
profiles as a sum of seeing convolved r 1 / 4 and exponential 
components. 

3.3. Isophotal analysis and photometric profiles 
NGC 4816: 

Fig. [|a shows the surface brightness profiles as well as the 
results of the isophotal analysis for the HST V band. We 
appended the "colour shifted" MDM R band data in the 
outer part of the galaxy (> 20"). The best fit of the MDM 
surface brightness profile is given by a sum of an r 1 / 4 law 
and an exponential component and is also shown in the 
figure. Inside the kinematically detected break radius (4" ), 
the surface brightness from the HST data shows the ex- 
ponential profile typical for a disk. Additionally the 04 
parameter is positive inside the core region. Outside the 
core the amplitudes of 04 parameter become slightly nega- 
tive. The ae parameter is slightly negative inside the core, 
while the amplitudes for the other Fourier coefficients are 
small and the mean ellipticity e = 0.2. We also show the 
surface brightness profile of the best model disk from the 
decomposition. For the remaining bulge we plot the 04 and 
parameter, the P.A. and e in Fig. ^b together with the 
data values for the whole galaxy. The corrected photomet- 
ric parameter for the central disk we determined from the 
HST data are presented in Table [|. Additionally mj- and 
r e derived from the MDM data for the whole galaxy and 
those available from the literature are listed in Table [| 
The discrepancy between the r e we measured and the lit- 
erature values is due to different maximal radial extension 
of the profiles. Our surface brightness profile extends to 
w 170", while the profiles of Andreon et al. (1996, 1997), 
Lucey et al. (1991) and Burstein et al. (1987) only reach 
to « 30" , 50" and 40" , respectively. Additionally only a 
pure r 1 / 4 profile was fitted, taking no account of a possible 
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NGC 4816: r =55.75" 



o 



m 



o 
o 



o 
o 



o 
o 



o 
o 



o 
o 




b) 



1.0 2.0 3.0 

a'/ 4 [arcsec 1 /*] 



1 1 1 1 I 



' "V bulge 

i I I I i I I I i I i i 



am mm - 



0.5 1.0 1.5 2.0 

a'/ 4 [arcsec'/ 4 ] 



2.5 



< 

CL 



I 



O 

o 



o 
o 



D 



o 
o 



< 

CL 



-Q 
I 



120 
100 
80 

60 
0.4 
0.2 

0.0 

2 

1 


-1 

-2 
2 
1 


-1 
-2 

2 

1 


-1 

-2 



140 
120 
100 
80 
60 

0.3 
0.2 
0.1 
0.0 
-0.1 



* 



I I 



III 



I I I I 



I I I 



LU dfe 



i i i i i i i i i i i i i 



i i i i i i i i i i i i i 



# * 



ill ill i i i i l i i 



1.0 2.0 3.0 

a'/ 4 [arcsec 1 ^] 



1 I 1 1 1 1 I 1 1 1 1 I 



I I I I I I I I I I I I I 



I I I I 



■ ■■■ 



0.5 1.0 1.5 2.0 2.5 

a'/ 4 [arcsec 1 / 4 ] 



Fig. 4. (a) Results of the isophote analysis of NGC 4816 from the HST V band data (open symbols) are plotted versus the 
major axis a. "Colour shifted" MDM R band data are appended for the outer part of the galaxy (a > 20"). The best fit of 
the MDM surface brightness profile - a sum of an r 1//4 law and an exponential component - is overplotted (fat solid line). The 
fit is only shown in the outer region, since the groundbased MDM profile starts to deviate from the HST profile inside 6" due 
to different PSFs. Left panels (top to bottom): (1) Surface brightness profile SB of the whole galaxy and the central disk only 
(dashed line), (2-4) Fourier coefficients 03/a, 04/a, ae/a. Right panels (top to bottom): (1) Position angle PA, (2) ellipticity e, 
(3-5) Fourier coefficients 63/a, 64/a be/a. (b) Results of the decomposition for the HST data (open symbols - whole galaxy; 
filled symbols - bulge only) plotted versus the major axis a. Left panels (top to bottom): (1-2) Fourier coefficients 04/a, ae/a. 
Right panels (top to bottom): (1) Position angle PA, (2) ellipticity e. 
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Galaxy 


Filter 




rh 


i 


PA d 


L d 


(D/B)tot 


D/B(r < r b ) 






[mag] 


["] 






[1O 9 L ] 






NGC 4816 


V 


17.79 


0.94 


64° 


100° 


1.26 


1.0 % 


6.9 % 


IC 4051 


V 


19.12 


0.74 


65° 


90° 


0.37 


0.4 % 


2.5 % 


IC 4051 


I 


17.53 


0.90 


65° 


90° 


0.79 


0.5 % 


2.8 % 



Table 2. Photometric parameters for the central disks of NGC 4816 and IC 4051 



exponential component. 



Galaxy 


Filter 


hit 
[mag] 


r e 
["] 


Source 


NGC 4816 


R 


11.61 


55.75 


this work 


NGC 4816 
NGC 4816 
NGC 4816 
NGC 4816 
NGC 4816 


R 

V 
V 
V 
B 


11.62 
12.66 
12.90 
12.82 
13.79 


27.54 
25.12 
17.14 
21.20 
21.20 


Andreon et al. (1996) 
Andreon et al. (1997) 
Lucey et al. (1991) 
Burstein et al.(1987) 
Burstein et al.(1987) 


IC 4051 


R 


12.46 


15.56 


Saglia et al. (1997b) 


IC 4051 
IC 4051 
IC 4051 
IC 4051 
IC 4051 
IC 4051 


R 

r 
V 
V 
B 
B 


12.77 
12.82 
13.10 
13.04 
13.92 
14.01 


10.72 
19.50 
18.66 
20.30 
22.39 
20.30 


Andreon et al. (1996) 
J0rgensen et al. (1992) 
Lucey et al. (1991) 

Burstein et al.(1987) 
J0rgensen et al. (1992) 

Burstein et al.(1987) 



Table 3. Photometric parameters from literature (and own 
photometry) for NGC 4816 and IC 4051 



IC 4051: 

The surface brightness profiles and the results of the 
isophotal analysis are shown in Fig. ||a for both, the John- 
son I and the V band exposure. We appended the "colour 
shifted" R band data from Saglia et al. (1997b) in the 
outer part of the galaxy (> 20"). The best fit of the R 
band surface brightness profile is given by a sum of an 
r 1 / 4 law and an exponentiell component and is also shown 
in the figure. In both the I and the V band the 04 parame- 
ter is positive inside 5" , which is the size of the co-rotating 
decoupled component, detected in the velocity curve. Ad- 
ditionally, the ag parameter is slightly negative inside the 
core, while the amplitudes for the other Fourier coefficients 
are small and the mean ellipticity e = 0.2. We also show 
the surface brightness profile of the best model disk from 
the decomposition. We plot the 0,4 and ag parameter, the 
P.A. and e for the remaining bulge in Fig. ||b together with 
the results for the whole galaxy. The corrected photomet- 
ric parameter for the central disk we determined from the 
HST data are presented in Table [2] for both filters. Ad- 
ditionally, the total magnitude m-r and the effective ra- 



dius r e for the whole galaxies available from the literature 
are listed Table ||. The central colour (V-I) is 1.355 mag 
and the V-I colour profile is show in Fig. |[ We derived 
a mean colour gradient A(V ~ I)/ Alog(r) within 10" of 
-0.018 mag/dex by a linear fit. For comparison, J0rgensen 
et al. (1992) found a B-r gradient of -0.09 mag per dex in 
radius inside r e . 

IC 4051 



L.36 



.34 



1.32 - 



> 



] .3 



' .38 



: .26 




1 



log (r/") 



Fig. 6. V - I colour profile of IC 4051 and the results of the 
least square fit for the colour gradient within 10" . 



4. Results 

The early - type galaxy NGC 4816 hosts a decoupled 
counter rotating core with a size of ib = 4" (Fig. ||a). In- 
side this core the H3 parameter, which measures the asym- 
metric deviations from a gaussian velocity distribution, is 
slightly positive where v < 0. This confirms Bender et al. 
(1994), who find that H 3 becomes systematically negative 
with increasing v/< a >. Hence, the LOSVD in the core 
of NGC 4816 are intrinsically asymmetric, showing an ex- 
cess of fast rotating stars. 

The isophotes inside 4" are clearly disky, which is derived 
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IC 4051; r e =15.56" 

a) 




a 1 /* [arcsec 1 / 4 ] a 1 / 4 [arcsec 1 / 4 ] 



Fig. 5. (a) Results of the isophote analysis of IC 4051 from the HST V and I band data (open symbols) are plotted versus the 
major axis a. "Colour shifted" R band data from Saglia et al. (1997b) (stars) are appended for the outer part of the galaxy (a > 
20"). The best fit of the R band surface brightness profile - a sum of an r 1//4 law and an exponential component - is overplotted 
(fat solid line). The fit is only shown in the outer region, since the groundbased R band profile starts to deviate from the HST 
profile inside 3" due to different PSFs. Left panels (top to bottom): (1) Surface brightness profile SB of the whole galaxy and the 
central disk only (short and long dashes), (2-4) Fourier coefficients as/a, 04/a, as/a.. Right panels (top to bottom): (1) Position 
angle PA, (2) ellipticity e, (3-5) Fourier coefficients 63 /a, &4/a be/ a., (b) Results of the decomposition for the HST data (open 
symbols - whole galaxy; filled symbols - bulge only) plotted versus the major axis a. Left panels (top to bottom): (1-2) Fourier 
coefficients 04/a, ae/a,. Right panels (top to bottom): (1) Position angle PA, (2) ellipticity e. 
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from the positive a4 parameter from the isophotal shape 
analysis (Fig. ^a). Since the disky isophotes extend to the 
same radius as the counter rotating core, detected in the 
velocity curves, it is evident that this component is a cen- 
tral stellar disk. Additionally, the decomposition of the 
surface brightness profiles into a central exponential disk 
profile and a bulge following an r 1 / 4 law, yields a perfect 
remaining bulge with 0,4 = (Fig. ^b). The central disk in 
NGC 4816 contributes only 1 % to the total V band light 
of the galaxy. This is even less than found for the kine- 
matically decoupled central disk inside NGC 5322 (Ben- 
der et al. 1994), which has D/G = 3 % (Scorza & Bender 
1995). Inside the break radius r& = 4" the central disk of 
NGC 4816 contributes 6.9 % to the total light. For a thin 
disk one would expect a positive a§ parameter from the 
isophote analysis (Scorza & Bender 1995), but not a neg- 
ative one as found in our data. Thus, with the assumption 
of having a thin central disk in NGC 4816 the "measured" 
inclination is only a lower limit. 

The slope of the Mgb line strength profile steepens inside 
the break radius r& (Fig. [2)3), which points to a metal en- 
riched core and agrees with results found for other pecu- 
liar cores by Bender & Surma (1992). The mean iron index 
<Fe> and the combined index [MgFe] do also increase in- 
side 8" . In the rotation curve the core is only showing up 
at 4" . As it is counter - rotating it has to compete against 
the superposed rotation of the bulge further out and be- 
comes only dominant inside 4" . In fact in (Fig. ^|a) you 
can see that the rotation curve seems to be influenced by 
the core out to 6" or 7" . Comparisons with stellar popru 
lation models (Worthey 1994) show that the metallicity □ 
of the core is Z = log ((Fe/H)/(Fe/H) ) » 0.25 dex super 
solar (Fig. 0a), which is comparable to values measured in 
central regions of normal galaxies (Faber et al 1995). On 
the other hand the metallicity of the hosting galaxy drops 
to values around solar metallicity (outer most data point 
at « 0.5 r e ). In the core as well as in the outer parts of 
the galaxy the dominating stellar population is old (8-12 
Gyr) (Fig. |a). 

From other studies it is already known that luminous, 
massive elliptical galaxies can reach a Mgf, to <Fe> over- 
abundance [Mg/Fe] « 0.4 dex super solar, while faint 
early - type galaxies have ks solar abundances (Gonza- 
les 1993; Fisher et al. 1995). NGC 4816 shows extreme 
Mgb to <Fe> overabundance of w 0.5 dex (Fig. 0b) with 
no significant abundance difference between the core and 
its surrounding, which is in agreement with other inves- 
tigations (Davies et al. 1993; Paquet, 1994). In a recent 
study, Thomas et al. (1997) show that in addition to a 
top heavy initial mass function (IMF), a short star for- 
mation time scale is needed to explain this light element 
enrichment in massive elliptical galaxies with an excess of 



4 Following Faber et al. (1995) we used the combined index 
[MgFe] = (Mgt <Fe>) 1//2 , which represents the mean metal- 
licity best. 



Supernovae II rate. 

For IC 4051 there is evidence for a co-rotating pecu- 
liar core with a radial size of rj, = 5" (Fig. |3]a) . In contrast 
to the core NGC 4816 it is not clear whether this central 
component in IC 4051 is really decoupled. From the pre- 
sented kinematic data (Fig. ^a), we cannot exclude that 
IC 4051 is simply a disky elliptical or SO galaxy being a 
superposition of a hot (rotating) bulge and a cold, fast ro- 
tating disk (Scorza & Bender 1995). In combination with 
the brightness profiles we derived from the HST images, 
we can compute the anisotropy parameter (Binney 1978, 
Kormcndy 1982): 

/ Vr V _ Vr/O'm ( ^ 

\a m ) ~ y/e/0.-e)' 

where v r is the rotation velocity, e the ellipticity 
and a rn the averaged velocity dispersion. Objects with 
(v/a)* > 0.7 are mainly rotationally flattened, while those 
having smaller values are flattened by anisotropic velocity 
dispersions (Davies et al. 1983). Since the anisotropy pa- 
rameter for IC 4051 is « 0.5 inside the "core" and w 0.4 
at r e it is rather unlikely that it is an SO galaxy. 
From the isophotal analysis and the photometric decom- 
position (Fig. |^) we inferred that this core is a central 
stellar disk. Like the disk in NGC 4816 it only contributes 
0.4 % to the total V band light of the galaxy and 0.5 % 
to the total I band light. Since the disks in "normal" SO 
galaxies contribute r* 4 - 11 % to the total V band light 
(Scorza & Bender 1995) IC 4051 is probably not an SO 
galaxy, but contains a decoupled central disk. Inside the 
break radius rj, = 5" the disk of IC 4051 contributes 2.4 % 
and 2.7 % in the V and I band, respectively. The V-I colour 
gradient inside the IC 405 l's core is shallow and in agree- 
ment with Carollo et al. (1997) it does not indicate any 
significant difference between the core and the surround- 
ing galaxy. 

In contrast to NGC 4816 the slope of the Mgh line profile 
does not increase inside the core of IC 4051 (r < 5" ) but is 
constant (Fig. ||b). This lack of steepening inside the core 
could be due to [NI] emission that can contaminate the Mg 
absorption feature (Goudfrooij & Emsellem 1996). On the 
other hand the gradient of the mean iron index <Fe> is 
enhanced towards the center starting at the break radius, 
pointing to a metal enriched core as found in NGC 4816. 
Comparisons with stellar population models (Fig. |8^) 
show that the metallicity drops from super solar inside 
the core (ss 0.25 dex) to -0.5 dex in the outer part of the 
galaxy (outer most data point at w 1 r e ). This metallicity 
difference of 0.75 dex between the central and the embed- 
ding galaxy is unusually high. The stars in the core are 
old (12 -17 Gyr), while an age estimate in the outer parts 
is not possible since the absorption lines are probably 
contaminated by emission like features in the outer parts 
of the galaxy. As in NGC 4816 the dominating stellar pop- 
ulation is highly [Mg/Fe] overabundant. 



Mehlert et al.: The kinematically peculiar cores of NGC 4816 and IC 4051 



11 
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Fig. 7. Measured absorption line indices of NGC 4816 plotted into the model grids of Worthey (1994); dashed lines - constant 
age in Gyr; dotted lines - constant metallicity log(Z/Zo). Filled symbols represent measurements inside the break radius r;,, 
open symbols those for measurements outside r b . (a) The age indicating H/3 versus the metallicity indicating combined index 
[MgFe]. (b) <Fe> versus Mg6 index. Obviously the models cannot reproduce the measured element ratios. 
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Fig. 8. Measured absorption line indices of IC 4051 plotted into the model grids of Worthey (1994); dashed lines - constant age 
in Gyr; dotted lines - constant metallicity log(Z/Z Q ). Filled symbols represent measurements inside the break radius rj,, open 
symbols those for measurements outside ib- (a) The age indicating H/3 versus the metallicity indicating combined index [MgFe]. 
Outside Tt the H/3 indices are emission contaminated and hence not reliable for age determination, (b) <Fe> versus Mgb index. 
Obviously the models cannot reproduce the measured element ratios. 
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5. Discussion 



open question, which can not be answered with our data. 



NGC 4816 contains a kinematically decoupled central 
disk, which is dominated by an old, metals enriched, highly 
[Mg/Fe] overabundant stellar population and is only con- 
tributing w 1 % to the total light of the galaxy. For the 
formation scenario of this disk we favour dissipational ma- 
jor merging for the following reasons: 

• A merging event is necessary to account for the kine- 
matic decoupling between core and bulge (Barnes 1992b). 

• Dissipational merging is required in contrast to dissipa- 
tionless processes because: (1) In a dissipational merger 
the gas can get to the center and form a flattened, disk - 
like core by dissipation and cooling (Barnes & Hernquist 
1996). Flattened systems like the core of NGC 4816 cannot 
be built through dissipationless processes. (2) While dis- 
sipationless merger tend to weaken metallicity gradients 
(White 1980), only dissipative processes can built up these 
gradients. Previously enriched gas settles in the core and 
may trigger the formation of metal enriched stars (Bender 
& Surma 1992) . The enhancement of the stellar metal- 
licity we find in the central disk of NGC 4816 require gas 
dissipation. 

• The accretion of low-mass gas rich galaxies can be ruled 
out since they are generally too metal poor (Bender & 
Surma 1992). To obtain the high core metallicity measured 
in NGC 4816 a major merging of massive, star dominated, 
but gas-rich pro-genitors is needed. 

• Finally, merging or interacting galaxies experience a 
burstlike star formation, that may produce flatter IMFs 
(Wright et al. 1988; Silk 1992) and hence could help to 
produce the [Mg/Fe] overabundance we measure in the 
stellar population of NGC 4816. 

Since the stellar population of the central disk is likely 
old, this merging events must have occured at high red- 
shift in the early phase of formation. However, the age 
indicating H^ line strength is strongly sensitive to H^ 
emission from younger populations. Hence, the estimated 
age is only an upper limit to the true age of the stars in 
the central disk. If the true age was significantly younger, 
the central disk in NGC 4816 could also have formed in 
late, gas-rich mergers, which are e.g. observed in ultra- 
luminous IRAS galaxies (Solomon et al. 1992, Prestwich 
et al. 1994). However, late mergers may have problems in 
producing the light element overabundances. 

The central disk in IC 4051 is also dominated by an 
old, metal rich, highly [Mg/Fe] overabundant stellar pop- 
ulation and contributes only « 1 % to the total light of 
the galaxy. In contrast to NGC 4816 the interpretation of 
the data with respect to the formation of the galaxy and 
its peculiar core is less evident. Whether the disk in IC 
4051 is really kinematically decoupled and hence whether 
this galaxy is indeed a major merger remnant remains an 
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Fig. 9. Distribution of the break radii r^ of all known kinemat- 
ically distinct cores. The positions of IC 4051 and NGC 4816 
are indicated by crosses. 



Table [| contains the 17 known early - type galaxies 
with clear evidence for a kinematically decoupled core. 
Most of the other cores have sizes smaller than 2 kpc 
(see the histogram Fig. ^) . With a break radius of r& = 
2.7 kpc NGC 4816 hosts one of the largest kinematically 
decoupled core ever detected in an early - type galaxy. 
Only the decoupled central of NGC 7626, the dominating 
galaxy of the Pegasus cluster (Jedrzejewski& Schechter 
1988; Balcells & Carter 1993) has a similar size. The 
peculiar core in IC 4051 (r^ = 3.4 kpc) is even bigger. 
To investigate whether the large cores of NGC 4816 and 
IC 4051 are systematically different to other, smaller kine- 
matically decouples cores, we compared the properties of 
their corresponding central disks. HST images are avail- 
able for 13 of the 15 galaxies that are not in Coma and 
are listed in Table ||[ We applied the D/B decomposition 
of Scorza & Bender (1995) to derive the properties of 
the central disks, that are present in 8 of these galax- 
ies. Since 4 of these disks were contaminated with dust 
we are left with reliable disk properties for NGC 1427, 
NGC 4365, NGC 5322 and NGC 5982. Plotting the cen- 
tral surface brightness SB versus the scale length r^ 
shows that the disks of NGC 4816 and IC 4051 simply rep- 
resent the "faint", "large" end of a sequence. We therefore 
speculate that the large cores of NGC 4816 and IC 4051 
had formation histories similar to the ones of the smaller, 
already known decoupled cores. Furthermore, the param- 
eters r/i and SB derived for the central disks of galaxies 
with cores are similar to the ones derived for "disky ellip- 
ticals" (Scorza & Bender 1995). 

Excluding Coma's two dominant galaxies (NGC 4874 & 
NGC 4889) and the central galaxy (NGC 4839) of the 
Coma substructure in the S/W (Colless & Dunn 1995), 
IC 4051 is the third brightest galaxy in the central square 
degree of the cluster (see J0rgensen & Franx 1994). 
NGC 4816 is even 0.75 mag brighter than IC 4051 in the 
R band and hence, both galaxies have luminosities above 
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Table 4. Galaxies with kinematically decoupled cores mentioned in the literature. We used H = 50 km / (s Mpc) to the 
determine the distances. At the "break radius" r& the decoupled core starts to dominate the rotation curve. Galaxies with an * 
belong to the Virgo cluster. 



L* . While NGC 4816 lies 49.0' S/W of the clusters center 
and only 21.8' apart from NGC 4839, IC 4051 lies only 
15.5' N/E of the clusters center. Whereas NGC 4816 is 
the largest host galaxy (r e = 37.8 kpc) of a distinct core, 
the size of IC 4051 is comparable the sizes of the host 
galaxie listed in Table ^. However, like the other galaxies 
with kinematically decoupled cores both galaxies do not 
deviate from the Fundamental Plane or the Mg - a rela- 
tion. Our total sample of 35 Coma galaxies is complete for 
galaxies brighter than Mg = -21.6 mag (14 galaxies) and is 
1/3 complete in the magnitude range -21.6 mag < Mb < 
-20.2 mag (21 galaxies). NGC 4816 and IC 4051 both be- 
long to the complete sample. Similarly, in the Virgo cluster 
there are only two galaxies (NGC 4406 and NGC 4472) 
that host a decoupled core and are brighter than -21.6 
mag. The core of NGC 4406 contributes w 20 % to the 
the V band light inside the kinematically break radius i'/, 
(Surma 1992), which is consistent with the central disks 
in the Coma galaxies. In the magnitude range -21.6 mag 

< Mb < -20.2 mag one galaxy with a decoupled core is 
known in the Virgo cluster (NGC 4365), but none found 
in the Coma cluster yet. In their detailed discussion of 
NGC 4365 's stellar population and kinematic Surma & 
Bender (1995) conclude that its core also must have built 
through dissipational major merging processes. This core 
also shows disk - like structure, its mean stellar population 
has a mctallicity which is sa 0.4 dex above the Sun and 
an age of « 7 Gyr. Inside the kinematically break radius 
r& it contributes « 20 % to the V band light. Note that 
the break radius r;, of the "Virgo cores" would compare to 

< 1.3" at Coma's distance. Thus with our observational 



setups (spatial resolution w 1.5") we missed all possibly 
existing cores with comparable size to the Virgo decou- 
pled cores. On the other hand there is no core detected in 
the Virgo cluster that has a similar size like those found in 
the Coma cluster galaxies NGC 4816 and IC 4051. Based 
on our data, we cannot derive any statistically significant 
difference between the Coma and the Virgo cluster with 
respect to the fraction of galaxies with kinematically de- 
coupled/peculiar cores. 
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